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r !p ■ Abstract 

We consider the quantum discord in a two-particle spin-1/2 system (dimer) gov- 
erned by the standard multiple quantum (MQ) Hamiltonian and reveal the relation 
\ between the discord and the intensity of the second-order MQ coherence. Since the 
coherence intensities may be measured in experiment, we obtain a method of de- 



tecting the quantum discord in experiment. Comparison of discord and concurrence 
is represented. 
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. 1 Introduction 



For a long time the quantum entanglement [Tf2"f3"ll4"ll5] has been considered as a proper measure 
of quantum correlations. However, it seemed out that entanglement does not cover all quantum 
correlations which are responsible for advantages of quantum processes in comparison with their 
classical counterparts. For instance, the quantum non-locality without entanglement is studied 
in [61T71I8] . the quantum speed-up with separable states is demonstrated in [911101111111211131114] . 

Quantum discord was first introduced in ref.[15] and then it was studied in refs. |16|fT7][T8] as 
an alternative measure of quantum correlations. This measure involves larger class of quantum 
correlation in comparison with entanglement. In particular, it may be nonzero even for the 
states with zero entanglement. Therefore it was noted [H] that namely discord might be re- 
sponsible for advantages of quantum information devices, in particular, for quantum speedup. 
The concept of quantum discord is intensively developing during last years. For instance, the 
operational interpretation of quantum discord in terms of state merging [19] and unitary invari- 
ant modification of discord [20] have been proposed. However, only very special cases have been 
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treated analytically [2T|22|23j . Moreover, relation between discord and physical characteristics 
of quantum system is not well studied yet. 

This paper is devoted to the investigation of discord in a system of two spin-1/2 particles 
(dimer) on the preparation period of the MQ NMR experiment [2^f25] with the thermodynamic 
equilibrium initial state in the strong external magnetic field. We show that the discord can 
be considered as a function of the second-order coherence intensity similar to the entanglement 
[26J. This conclusion gives us a method to measure the discord in the MQ NMR experiment. 
Thereby one should note ref.[27] where the magnetic susceptibility has been introduced as 
a physical characteristic of a quantum system allowing one to measure the entanglement in 
nitrosyl iron complexes, as well as ref. [28] where the relations between discord and such 
physical characteristics as internal energy, specific heat, and magnetic susceptibility have been 
studied. 

We derive the explicit formula for the discord in a dimer with the dipole-dipole interaction 
(DDI) governed by the standard MQ NMR Hamiltonian [25] in SecfSJ The dependence of the 
discord on the second-order MQ NMR coherence intensity is studied in Secj3) The comparison 
with analogous dependence of the concurrence (as a measure of entanglement) is given in the 
same section. Conclusions are given in SecJH 



2 Discord in the system of two spin-1/2 particles with MQ Hamiltonian 



The dynamics of a system of two spin-1/2 particles with the DDI on the preparation period 
of the MQ NMR experiment [24125] in the strong external magnetic field is governed by the 
standard MQ NMR Hamiltonian Hmq, which reads 

HMQ = j(lt4 + Ill2) (1) 

in the rotating reference frame [29], where D is the coupling constant between two spins, which 
reads D = — (1 — 3 cos 2 6*12) in the case of DDI. Here 7 is the gyromagnetic ratio, ri2 is the 

r 12 

distance between spins, 812 is the angle between the vector r*i2 and the external magnetic field 
H ; If = Ij X + iljy and I J = Ij x — ilj y (j = 1,2) are the raising and lowering operators of 
spin j; I a j are operators of the jth spin angular momentum projection on the axis a — x,y, z. 
Starting with the thermodynamic equilibrium initial state of the dimer in the strong external 
magnetic field, 



( where I z = l\ z + I<i z is the ^-projection of the total spin angular momentum , /3 = hu>o/kT 
is the dimensionless inverse temperature, uq = 7|i7o|) we obtain the following evolution of the 
density matrix [26] (we use the standard basis of the eigenstates of the operator I z : 1 00) , |01), 
|10), |11)): 
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which is the solution to the Liouville equation (K — 1) 
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= I h mq,P(t)]. 



Note, that the density matrix ([31) has the structure 
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which is a particular case of the more general X-matrix whose discord has been analytically 
studied in [22]. In our case, formulas derived in [22] may be simplified yielding explicit ex- 
pression for the discord Q as a function of (3 and r. By definition [T7|16f22] the discord is a 
difference of two components: the total mutual information I(p) and the classical correlations 
C(p), i.e. 

Q(p) = I(p)-C(p). (6) 
Here the total mutual information I(p) is given by [17] 



J(p) = S(p A ) + S(p B ) + £ A, log 2 A„ 



(7) 



i=o 



where A and B mark the first and the second spins respectively, p A = Tr^p = diag(pn + 
P22,P33 + P44), P B = TtaP = diag(pn + p 33 , p 22 + P44) are the reduced density matrices and A.,- 
are the eigenvalues of the density matrix (111): 



^ _ cosh (3 + sinh /3 , _ cosh /3 — sinh (3 
~ 2(1 + cosh (3) ' 1 ~ 2(1 + cosh (3) 



Xo — X: 



2(1 + cosh P) 



Quantities S(p A ) and S(p B ) are the appropriate entropies [TTJ (S(p) = — Tr(plog 2 p)) which 
read in our case 
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S(P A ) = S(p B ) = (9) 

1 (cosh/3 + l) 2 - £ 2 sinh 2 /3 £sinh/3 cosh f3 + 1 + f sinh (3 

~2 ° g2 4(1 + cosh (3) 2 2(1 + cosh J) ° g2 cosh + 1 - £ sinh /3 

where we introduce the time-dependent parameter £ = | cos(.Dt)|, < £ < 1. Following ref . |22] . 
we assume that the projective measurements are performed over the subsystem B. In our case, 
the expression for the classical correlations C(p) (see refs. [TTf2"2"] ) in Eq.Q is simplified: 

C(p) = S(p A )- min flfo, AO, (10) 

r?={0,l} 

where 

nfa,/U)=PoS'o+PiS'i (11) 

with 



$ = 2 — g2 — 2 2 — g2 — 2 — ' ^ ' 

Pi = l( 1 + (-l)^(2(pu + Pas) - 1)) , (13) 

^ = i[(l-r, 2 )|p 14 | 2 + (14) 

^ (2(pn + p 22 ) - 1 + (-1)^(1 - 2(p 22 + p 33 ))) 'J 1 \ z = 0, 1. 



Here 77 is an arbitrary parameter, < rj < 1. The minimum of /3, r) in eq.f lTO]) corresponds 

tO T) = 0, 



^0,/3,0 = log 2 (l + e^)- ln2( / j^ e/3) , (15) 

which is proven in Appendix. Thus, the optimization problem, which is inherent in the definition 
of the quantum discord [IT], is completely solved for the dimer in the MQ NMR experiment. 
Remember, that optimization problem for general X-matrix [22] is reduced to the minimum of 
six values. 

Now we are able to write a single expression for the discord. In fact, combining eqs. (jnHHII]) 
and taking into account that Yh=o \ l°g 2 \ = ~ 2 f2(0, j3, £) (which can be simply derived sub- 
stituting Aj from Eqs.flU) into the expression Y%=o A, log 2 Aj) we obtain the following expression 
for the quantum discord: 



Q(p,o = s( P A )-m^o= (is) 
io g2 (i + e?) - ^rr+j) ~ \ log2 ( (cosh + 1)2 _ e sinh2 

£ sinh (3 cosh /3 + 1 + £ sinh /3 

log 2 
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Fig. 1. The discord Q as a function of dimensionless inverse temperature f3 and time-dependent pa- 
rameter £ 

The three-dimensional plot of the discord Q (fl6l) as a function of dimensionless inverse tem- 
perature j3 and time- dependent parameter £ is depicted in FigJT] 



3 Quantum discord and concurrence as functions of second-order coherence in- 
tensity 

Since the discord has been introduced as a measure of quantum information, it is important to 
have a tool allowing one to measure the discord in experiments [28]. Regarding the MQ NMR 
experiment, the MQ coherence intensities should be taken as measurable quantities for this 
purpose [2lf25] . 

Remember, that the standard MQ NMR experiment consists of four periods [25]: (i) the prepara- 
tion period, where the Hmq Hamiltonian is generated by the special sequence of radio-frequency 
pulses; duration of this period is r, (ii) the evolution period, where the system evolves with- 
out radio-frequency pulses; the effective Hamiltonian on this period is AI Z , where A is some 
constant parameter; duration of this period is t, (hi) the mixing period, where the system is 
irradiated by the special pulse sequence generating the —Hmq Hamiltonian; duration of this 
period is r, and (iv) detection. Due to the special structure of the MQ Hamiltonian Hmq, the 
transitions between eigenstates with different ^-projections of the total spin angular momen- 
tum become possible. The longitudinal magnetization after the mixing period involves signal 
generated by these transitions. The MQ NMR experiment allows one to split the longitudinal 
magnetization into so-called MQ NMR coherences of different orders k. The /cth-order MQ NMR 



5 



coherence combines contributions from the transitions between states whose z-projections of 
the total spin angular momentum differ by k. Such transitions are possible in clusters of at least 
k correlated spins which are formed due to the DDI interactions in the MQ NMR experiment. 
We do not represent details of the MQ NMR experiment, which may be found, for instance, 
in [23] • We only remind how different MQ coherence intensities G k (r) may be selected from 
the Fourier transform of the longitudinal magnetization (I z ) after the mixing period of the MQ 
NMR experiment. 

The longitudinal magnetization (I z ) after the mixing period reads [29] : 

(J,(r,t))=Tr(p(r,t)J«), (17) 
where p(r, t) is the density matrix after the mixing period, 

p( T , t) = jHMQT e -iAI z t e -iH MQ T poe iH MQ r e iAI z t e -iH MQ T _ ^ 

Then Eq. (jl7p may be written as follows: 



N 



(I z (r,t)) =Tr {e~ iA ^p(T)e lAht p h \r)) = E e- ikAt G k (r), (19) 

k=-N 

where p is defined in Eq. (|3]) , while the density matrix p ht reads 



p ht ( T ) = e -iHMQTj ze iH MQ T^ ( 20 ) 

and G k is the fcth-order MQ coherence intensity, 

G k = p k p h \. (21) 
In Eq.( )2T]) . we use the following representations of the density matrices: 



N N 

P= E Pk, P M = E Pk, (22) 

k=-N k=-N 

where pk and p^ are those parts of the matrices p and p ht respectively which are responsible 
for the /cth-order coherence [33], i.e they satisfy the conditions 

e~ iAht p k e lAht = e~ lkAt p k , e^'Vf e iAJ '* = e~ lkAt p h k l . (23) 

Let us remember [3T|32|33] that only the zeroth- and ±2nd-order MQ NMR coherence inten- 
sities are generated in the case of dimer (which might be considered as a particular model 
with the nearest neighbor interaction) and the intensities G±2 for the state ([3]) are defined as 
[33135125] 
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G = G± 2 = Pu(t) P £(t) = ^tanh^sin 2 (Dr) = ^anh^(l - £ 2 ), (24) 

Now we proceed to study the relation between the discord Q and the coherence intensity G. 
We follow ref. [26], where the concurrence (as a measure of quantum entanglement) has been 
derived as a function of (3 and r, 



C^,0 = max(o, Vr ^^l "J' e=|cos(Dr)|, (25) 

and its dependence on the second-order coherence intensity G has been studied for dimer. We 
show that the second-order coherence intensity is an appropriate observable quantity, which 
may serve to detect not only the concurrence, but also the quantum discord in the MQ NMR 
experiment with dimers. 

We describe the relation between Q and G in two ways resulting to the discord as a function 
of either (3 and G (Sec J3.1j) or G and £ (Sec l3.2l) . We also compare these relations with the 
analogous relations for concurrence. 



3.1 Discord Q and concurrence C as functions of dimensionless inverse temperature (3 and 
coherence intensity G 



Using Eq. (1241) we may express £ in terms of G, 



« = f^Smr (26) 

and substitute it into eq. flTE]) to end up with the discord Q as the function of (3 and G: Q(/3, G). 
The explicit expression for Q(G,£) is very complicated and it is not represented here. Sub- 
stituting Eq. (l26p for £ into Eq. (1251) we obtain the concurrence as the function of [3 and G 



C((3, G) = max ( 0, y 2G tanh ^ - ^^r| | , (27) 

The concurrence C and the discord Q as functions of G at three different values of (3 are 
compared in Figj2^,. In this case, in accordance with Eq.f l26|) . G varies in the range 

< G < G« ,09) = ~ tanh(/?/2), (28) 



so that the maximal values of the concurrence and discord are C(/3,G$ ax ) and Q(/3,G^ LXJ 
respectively. 
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According to Eq. fl27|) . the concurrence is positive (i.e. the state is entangled) if (26] 



g£LG9) < G(r, (3) < G^lM, G^M = — — J— (29) 



unlike the discord Q, which is positive over the whole interval ( 1281) . 

The comparison of concurrence C and discord Q as functions of (3 at three different values of 
G is depicted in FigFSb. Here parameter j3 is restricted as 



Pinl(G)<f3, /3i 1 L(G) = 2tanh^ 1 (2G), (30) 

which follows from Eq. (l26l . Thus the minimal value of discord is Q(/3^J n ,G). Again, the con- 
currence is positive if 



max(/3i\ ) n (G'),/3il(G))</3, (31) 

where p^ in = 2tanh~ 1 (X 2 ) and X is the unique (for positive G) positive solution to the 
equation 



— + V2G X - - = 0, (32) 
2 2' v ' 



which is obtained from Eq.(j27l). Thus, the minimal value of the concurrence is C(max(Pl n - n , filnin)-, G) 
which is zero if < ^] n . 



Emphasize that the discord may be valuable even if the concurrence is zero, which is confirmed 
by Figs. Hi and [2b (G = 0.1). 



3.2 Discord Q and concurrence C as functions of the coherence intensity G and time- 
dependent parameter £ 



Similarly, we may express the dimensionless inverse temperature /3 in terms of the coherence 
intensity G from Eq. (l24|) . 



(3 = 2tanh- 1 (2G/(l -f 2 )), (33) 

and substitute it into eq. ljTBl) obtaining the discord Q as a function of G and £: Q(G,^). The 
explicit expression for Q(G, £) is not represented here, because it is too cumbersome. 

Eliminating (3 from Eq. (l25|) using Eq. (l33|) we obtain the concurrence as a function of G and £: 
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Fig. 2. The discord Q and the concurrence C as functions of the second-order coherence intensity 
G and inverse temperature /3. (a) The discord Q and the concurrence C versus the second-order 
coherence intensity G at different /3, /3 = 1,2,5. The appropriate maximal values of G are fol- 
lowing: G^L=0.231, 0.381, 0.493. At these G, discord and concurrence reach their maximal val- 



ues: Q r 



tive at G^] n =0.167, 0.029, 9.0 x 10 5 . (6) The discord Q and the concurrence C versus the in- 
verse temperature /3, at different G, G = 0.1,0.25,0.4. The appropriate minimal values of f3 are 
followin 
ues: Q r 



0.160, 0.473, 0.942, G n 



0.069, 0.552, 0.973. The concurrence becomes posi- 



following: /3^] n =0.405, 1.099, 2.197. At these j3, discord and concurrence reach their minimal val- 



At = 1-295 (for G = 0.1) and at /3^ n 



C(G,0 = max 



0.029, 0.189, 0.531, C max = 0, 0.125,0.620. The concurrence becomes positive at 

1.099, 2.197 (for G = 0.25, 0.4) 

2 J ' 



2\2 



(34) 



The discord and the concurrence as functions of G at three different values of £ are compared 
in FigJSK In this case G varies in the range 



(35) 



so that the maximal values of discord and concurrence are Q{G^ X ^) and C(G$ ax ,£). The 
almost linear behavior of the concurrence C is explained by Eq. (]34p . which means that C ~ G 
for G < 1. 

According to Eq. ([3~4~j) . the concurrence is positive (i.e. the state is entangled) if 



g2.«) < g < g2.(0, «SL(0 = j (d - - e - a - e?» 

unlike the discord, which is positive over the whole interval ( 135]) . 



(36) 



Comparison of the discord and the concurrence as functions of £ at three different values of G 
is depicted in FigJHJ). Here we have (according to eq. (!33~]) ) 



(37) 
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Entanglement witness, G f=|cos(2Ar)| 



Fig. 3. The discord Q and concurrence C as functions of the second-order coherence intensity and 
time-dependent parameter £. (a) The discord and concurrence versus G at different f = 0.9, y/2/2, 0. 
The appropriate maximal values of G are following: (2^=0.095, 0.250, 0.500. At these G, discord 
and concurrence reach their maximal values: Qmax — 

0.286, 0.601, 1, C max = 0.436, 0.707, 1. The 

(2) 

concurrence becomes positive at G min = 0.062, 0.129, 0.207. (6) The discord and concurrence versus £ 

(2) 

at different G = 0.1, 0.25, 0.4. The appropriate maximal values of £ are following: Q n ' ax =Q.S9A, 0.707, 
0.447. At these £, discord and concurrence reach their maximal values: Qmax — 

0.298, 0.601, 0.850, 

C max = 0.447, 0.707, 0.894. The concurrence becomes positive at f^L = 0.806 (for G = 0.1) and at 
£ = (for G = 0.25, 0.4) 

Thus the maximal values of the discord and concurrence are Q(G, anc ^ C(G, CmD respec- 
tively. The concurrence is positive if 



dl(G)<e<d 2 L(G), (38) 
where £mi n (G) = \/T— and X is the unique (for positive G) positive solution to the equation 



X 4 

— - 2GX 3 - 2G 2 = 0, (39) 
which is derived from Eq. f!34p . 

Similar to Figj2j Figs. [3^ and|3]3 (G = 0.1) demonstrate that the discord may be valuable even 
if the concurrence is zero. 

Thus, having functions Q(/3,G) and C(j3,G) (or Q(G, £) and C(G,£) ) we are able to find the 
discord and concurrence measuring the second order coherence intensity at the given inverse 
temperature /3 (or at the given value of the parameter £ = | cos(Dr)\). 

The graphs in Figs(H [2] and [3] demonstrate that the maximal values of both discord and con- 
currence equal to unity, which corresponds to £ = 0, G — 1/2 (or (3 — > oo, G = 1/2). 

Finally we remark, that the limit (3 — » oo (T = 0) yields the pure state. In fact, the density 
matrix ([3]) reads in this limit: 
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P(r) 



2 



(40) 















V 



ism(Dr) 1 




so that p 2 = p, which is the criterion of pure state. It is known [T3] that the discord coincides 
with the entanglement for this state. 



4 Conclusions 

We show that the quantum discord in dimers (whose spin dynamics is governed by the Hamil- 
tonian Hmq ([I]) ) may be expressed in terms of the second-order coherence intensity, similar to 
the entanglement [25] . Due to this fact, the discord, which is the function of the dimensionless 
inverse temperature and time (i.e. /3 and £ = | cos(-Dr)|), may be considered either as a func- 
tion of the dimensionless inverse temperature (3 and the coherence intensity G (i.e. Q((3,G)) 
or as a function of the second-order coherence intensity G and the time r (i.e. Q(G,£)). Thus, 
considering the MQ NMR experiment with dimer, one can detect the quantum discord in the 
sample measuring the second-order coherence intensity therein. 

Comparing the discord and the concurrence in FigsfS] and [3J we observe that there is some 
quantitative difference between discord and concurrence, while the general behavior of these 
two measures of quantum correlations is very similar. To the most important difference one 
must refer the fact that the concurrence may be zero for positive G, while discord vanishes 
only at G — 0. Even if concurrence is zero, the discord may be valuable. This confirms the 
advantage of discord as a measure of quantum correlations in a system. The fact that graphs of 
concurrence are above the graphs of discord in many cases (see FigsfSJ [3J) indicates that these 
measures take into account quantum correlations in different manners. 

One has to emphasize that the requirement to have non-zero discord is much less rigorous than 
the requirement to have non-zero entanglement in a quantum system, which is confirmed by 
Figs. [2] and [3j This enriches the variety of quantum systems applicable in quantum devises. 
In particular, one has to note that the non-zero discord has been observed in the liquid-state 
NMR J3HEZEB] , where the quantum entanglement is absent. This creates a new material basis 
for development of quantum information devises. 

Authors thank Professor E.B.Fel'dman for useful discussions. This work is supported by the 
Program of the Presidium of RAS No. 8 "Development of methods of obtaining chemical com- 
pounds and creation of new materials" . 



5 Appendix 

Let us show that 
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min (fi(0, 0, 0,0(1, ft £)) =«(0,/3,O- 



(41) 



The explicit expression for f2(0,/3,£) is given in Eq. (|T5|) . while f2(l,/3,£) reads: 



ft = i log 2 ((l + cosh 0) 2 - e sinh 2 0) - (42) 



cosh/3 



2(1 + cosh /3) 



log 2 (cosh 2 /3 -£ 2 sinh'/3) - 



£ sinh /3 ^ (1 + cosh /3 — £ sinh 0) (cosh /3 + £ sinh /?) 



2(cosh 0+1) 2 (1 + cosh /3 + £ sinh 0) (cosh /3 - £ sinh 0) ' 

We see that f2(0, /3, £) does not depend on £ and f2(l, /?, 1) = £7(0, /?, 1). Thus, in order to prove 
equality ( l4"T|) . it is enough to show that f2(l,/3,£) is the decreasing function on the interval 
< £ < 1 and consequently reaches its minimal value at £ = 1. 

For this purpose we calculate the derivative of eq. (l4*2|) with respect to £ obtaining the following 
result: 



^ fl(l 8 £) sinh/3 (1 + cosh 0) cosh — £ 2 sinh 2 — £ sinh 

~dX ^ ' ^ ^ ~ 2(cosh + 1) ° g2 (1 + cosh 0) cosh - £ 2 sinh 2 /3 + £ sinh /3 ' 



Simple algebraic estimations show that 

(1 + cosh /3) cosh /3 - £ 2 sinh 2 /? - £ sinh < ^ 
< ( 1 + cosh /3) cosh /3 - £ 2 sinh 2 /3 + £ sinh ^ ^ ^ 

on the interval < £ < 1. Therefore the derivative (1431) is negative. Consequently, function ( |42l) 
decreases on this interval and achieves the minimal value at £ = 1. Hence Eq. (j4Tj) is valid. 
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